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Abstract 

Designing mobiles to harvest ambient energy such as kinetic activities or electromagnetic radiation 
will enable wireless networks to be self sustaining besides alleviating global warming. In this paper, 
the spatial throughput of a mobile ad hoc network powered by energy harvesting is analyzed using a 
stochastic -geometry model. In this model, transmitters are distributed as a Poisson point process and 
energy arrives at each transmitter randomly with a uniform average rate called the energy arrival rate; 
upon harvesting sufficient energy, each transmitter transmits with fixed power to an intended receiver 
under an outage-probability constraint for a target signal-to-interference-and-noise ratio. It is assumed 
that transmitters store energy in batteries with infinite capacity. By applying the random-walk theory, the 
probability that a transmitter transmits, called the transmission probability, is proved to be equal to one if 
the energy-arrival rate exceeds transmission power or otherwise is equal to their ratio. This result and tools 
from stochastic geometry are applied to maximize the network throughput for a given energy-arrival rate 
by optimizing transmission power. The maximum network throughput is shown to be proportional to the 
optimal transmission probability, which is equal to one if the transmitter density is below a derived function 
of the energy-arrival rate or otherwise is smaller than one and solves a given polynomial equation. Last, 
the limits of the maximum network throughput are obtained for the extreme cases of high energy-arrival 
rates and dense networks. 

Index Terms 

Energy harvesting, mobile ad hoc networks, throughput, power control, stochastic processes, mobile 
communication 

I. Introduction 

Recent years have seen increasing popularity of mobile devices such as smart phones and tablet comput- 
ers, giving rise to two design issues among others. First, the power consumption of mobile-device networks 
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makes an escalating contribution to global warming. Second, conventional batteries that power mobile 
devices periodically interrupt their operation due to finite battery lives; battery recharging using power 
cords causes inconvenience to users. These issues provide strong motivation for powering mobile devices 
by harvesting ambient energy such as solar energy, vibration, kinetic activities and electromagnetic radiation 
[1]. The capacity of mobile-device networks powered by energy harvesting remains largely unknown, which 
is addressed in this paper. 

This paper considers a mobile ad hoc network (MANET) where transmitters are modeled as a homoge- 
neous Poisson point process (PPP). Energy arrives randomly at a transmitter with a fixed average rate, called 
the energy -arrival rate. The energy-arrival process is modeled as an independent and identically distributed 
(i.i.d.) sequence of random variables and different processes are assumed independent. Each transmitter 
deploys an energy harvester that stores arriving energy in a rechargeable battery. Upon harvesting sufficient 
energy, a transmitter transmits with fixed power to an intended receiver under an outage-probability 
constraint for a target signal-to-interference-and-noise ratio (SINR). Based on the above model, the network 
spatial throughput is maximized by optimizing transmission power for a given energy-arrival rate. 

A. Prior Work and Motivation 

The fluctuation in harvested energy due to random energy arrivals requires redesigning existing transmis- 
sion algorithms for wireless communication systems. For single-user systems with infinitely backlogged 
data, knowledge of channel states and the temporal profile of energy arrivals can be exploited for controlling 
transmission power to maximize the throughput [2], [3]. The optimal power-control policies are shown to 
be variations of the classic water-filling policy such that the causality of energy arrivals and finite battery 
capacity are accounted for. Single-user systems with both bursty data and energy arrivals are studied in 
[4] and the optimal policies for joint queue-and-power control are characterized. Adaptive transmission for 
broadcast channels with energy harvesting has been also investigated [5], [6]. In [5], the optimal power- 
control for a two-user single-antenna broadcast channel is shown to attempt to allocate a fix amount of 
harvested energy to the user with the better channel before giving the remaining energy to the other user. A 
two-user multiple-input-multiple-output broadcast channel is considered in [6] where one user receives data 
and the other scavenges transmission energy, and the precoder at the base station is designed to optimize 
the tradeoff between the data rate and the rate of harvested energy. 

Wireless networks with energy harvesting have been studied [7]-[9]. For a wireless sensor network 
with energy harvesting and based on a simple channel model that omits channel noise and path loss, the 
probability that a sensor successfully transmits a data packet to a fusion center is analyzed in [7] for 
different multiple-access protocols including time-division multiple access and Aloha like random access. 
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Managing traffic load in time and space is important for wireless networks to be self sustaining through 
energy harvesting. Targeting wireless sensor networks with energy harvesting, distributive strategies are 
proposed in [9] for adapting traffic load to the spatial-and-temporal energy profile and evaluated using 
a network prototype. For a two-user interference network with energy harvesting, the data and energy 
arrivals are modeled as Bernoulli processes and the stability region is characterized such that it comprises 
all data-rate pairs under the constraint of finite data-queue lengths [8]. In view of prior work, there are 
few results that quantify the tradeoff between the network throughput and the energy-arrival rate though 
such results specify the fundamental limit of the network performance. The aforementioned tradeoff is 
investigated in the sequel using a stochastic-geometry approach. 

Stochastic geometry provides a set of powerful mathematical tools for modeling and designing wireless 
networks [10]. MANETs based on random access and carrier-sensing multiple access have been modeled 
using the PPPs [11], [12] and Matern hard-core processes [13], respectively. Cellular networks have been 
shown to be suitably modeled using the Poisson Voronoi tessellation [14]. Models of coexisting networks 
can be constructed by superimposing multiple point processes [15], [16]. Stochastic geometry models of 
wireless networks have been deployed to quantify the network-performance gains due to physical-layer 
techniques such as opportunistic transmission [17], bandwidth partitioning [18], successive interference 
cancellation [19], and multi-antenna techniques [20]-[24]. The performance metric typically considered in 
the literature is the network spatial throughput under a constraint on the outage probability for a target 
SINR, which is also adopted in this paper. Using this metric, most prior work focuses on deriving the 
outage probability using techniques such as the Laplace transform [11], [14] and probabilistic inequalities 
[12], [17]. This paper considers a MANET with Poisson distributed transmitters similar to the existing 
literature (see e.g., [12]). However, the transmitters in the current network model are powered by energy 
harvesting instead of reliable power supplies as in prior work. The consideration of energy harvesting 
introduces several new design issues including the relation between the network throughput and energy- 
arrival rate, the corresponding optimization of transmission power, and the effect of finite energy storage, 
which are investigated in the sequel. 

B. Contributions and Organization 

For exposition, a few definitions and notations are provided as follows. Time is slotted. Define the 
transmission probability p as the probability that a transmitter transmits and the network interference 
temperature as the maximum active transmitter density under the outage-probability constraint. Let Ao 
denote the transmitter density, A e the energy- arrival rate, and Z the amount of energy harvested by a 
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typical harvester in an arbitrary slot. 

The main contributions of this paper are summarized as follows. 

1) Assume infinite battery capacity. Using the law of large numbers and random-walk theory, it is 
proved that p is equal to one if P is larger than A e or otherwise p is equal to X e /P. 

2) Consider the case of finite battery capacity. Bounds on p are derived, which converge to the results 
stated above as the battery capacity increases. Moreover, two special cases are considered. If Z is 
bounded and no larger than P, it is shown that p is equal to X e /P so long as the battery capacity 
is larger than IP. If Z is a discrete random variable, p is analyzed using Markov chain theory. 

3) Assume infinite battery capacity. By applying derived results on transmission probability and tools 
from stochastic geometry, the network throughput is maximized by optimizing P for given A e . 
Consider the condition that Ao is smaller than the network interference temperature evaluated for 
equal P and A e . If this condition holds, the maximum throughput R* is shown to be 

R* = A log 2 (l + c?) (1) 

where 6 is the target SINR. If the aforementioned condition is not satisfied, 

fl* = ^l°g 2 (l + *) (2) 
where the optimal transmission power P* is larger than A e and solves a derived polynomial equation. 

4) Furthermore, the limits of the maximum network throughput are obtained for the extreme cases of 
high energy-arrival rates (A e — > oo) and dense networks (Ao — > oo). Specifically. 

lim R*{X e ) = min ( A , ^) log 2 (l + 9) 

lim ir(A o ) = -^log 2 (l + 0) 
where p t is a positive constant determined by the maximum outage probability. 
The remainder of this paper is organized as follows. The network model and performance metric are 
described in Section II. The transmission probability is analyzed in Section III. The results are applied to 
maximize the network throughput in Section IV. Numerical results are presented in Section V followed 
by concluding remarks in Section VI. 

II. Model and Metric 

A. Network Model 

As illustrated in Fig. 1, the transmitters {T} of the MANET are distributed in the Euclidean plane R 2 
following a homogeneous PPP $ with density Ao, where T denotes the coordinates of a transmitter. Each 

'A typical point is selected from a spatial point process by uniform sampling. 
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Fig. 1. Single-antenna transmitters of the MANET are modeled as a homogeneous PPP in the horizontal plane. Each transmitter 
is powered by an energy harvester and transmits to an intended receiver at an unit distance. 

transmitter is associated with an intended receiver located at a unit distance. The signal transmitted by T 
with power P is received by a receiver located at X with power equal to P\X — T\~ a where a > 2 is 
the path-loss exponent. In other words, propagation is characterized by path loss while fading is omitted 
to simplify notation. 2 

Time is partitioned into slots of unit duration with t denoting the slot index. The amount of energy 
harvested by the typical harvester in the t-th slot is represented by the nonnegative random variable Z t . 

Assumption 1. The energy-arrival process {Z t } C II 4 " is an i.i.d. sequence and independent with other 
energy-arrival processes. Moreover, the cumulant generating function of the random variable (Z t — j3), 
namely In E [e^ Zt '^] , has a root r*(/3) > if > X e . 

This assumption allows the use of results on the large deviation of random walks in the subsequent 
analysis [25]. Let B denote the battery capacity identical for all harvesters. Since the typical transmitter, 
represented by To, transmits with fixed power P, the battery level of the corresponding (typical) harvester, 
denoted as St, evolves as 

St = mm(S t - 1 + Z t -PI{S t - 1 >P),B), t = l,2,--- (3) 

where Sq = and the indicator function 1(A) for an event A is equal to one if A occurs or else is zero. 

2 The consideration of fading has no effect on the analytical results except that the parameter /i e defined in (35) has to be 
redefined by including additional random variables modeling channel fading coefficients. 
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B. Performance Metric 

Assume infinitely backlogged and packetized data. A transmitter transmits one data packet whenever 
the corresponding battery level exceeds P. Thus, the transmission probability p can be written as 

1 n 

p = lim -VE[/(S t >P)]. (4) 

n=oo n L — » 
t=l 

According to Coloring Theorem [26], the process of active transmitters, denoted as 11, is a PPP with 
density At = p\o. Data is encoded at a fixed rate log 2 (l + 6) bit/s/Hz with 9 > being the target SINR. 
Correct decoding of a data packet requires the received SINR to be no smaller than 9 or else an outage 
event occurs. The outage probability P out is defined as P out = Pr(SINR < 9) where SINR represents the 
received SINR at the receiver for To. It is assumed that the receiver for To is located at the origin, which 
does not compromise the generality based on Slyvnyak's Theorem [27], and that noise has unit variance. 
Based on these assumptions, P out can be written as 

Po U t = Pr( Yl P\T\- a >j~l 
\ren\{T } 

\Ten / 

where (5) uses Slyvnyak's Theorem. To ensure the quality-of-service, an outage-probability constraint 
is applied such that P out < e with < £ « 1. The performance metric is the network throughput R 
(bit/s/Hz/unit-area) defined as 

i? = A t log 2 (l + 0) (6) 
= X plog 2 (l + 9) (7) 

where p is controlled by adjusting P such that the outage-probability constraint is satisfied. To be precise, 
R should be scaled by the success probability (1 — e) but this factor is a constant and thus omitted for 
ease of notation. 

III. Transmission Probability 

A. Infinite Battery Capacity 

Deriving transmission probability requires analyzing the distribution of battery levels at energy harvesters. 
By substituting B — > oo into (3), the battery level at the typical energy harvester with infinite battery 
capacity evolves as 

St = S t -i + Z t -PI{S t - 1 >P). (8) 
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The distribution of St can be related to the threshold-crossing probability for a random walk as follows. 
Denote the instants when the battery level crosses the threshold P from below as h,t2, - ■ ■ , namely that 
<St n -i < P and St n > P for n = 1, 2, • • • . These time instants are grouped into the set T = {ti, t2, ■ ■ ■ }. 
Moreover, define the random variable Z t = Z t — P and two random processes {Gt} and {G' t } as 

G t = m a x(G t -i + Z t ,0) (9) 
( S t , t€T 

G' t = { (10) 
(G' tl , tiT 

with Go = and G' Q = P. Based on (9), the probability that {Gt} crosses a threshold x > in the t-th 
slot can be written as 



Pv(G t > x) = Pr ( max ( Z t , Z t + Z t _i, ■ ■ ■ , ^ Z n 

V V n=i / / 

that is equivalent to the probability that the random walk {X)n=i ever crosses tne threshold x by the 

i-th step. Applying Kingman bound on the threshold-crossing probability for a random walk [25] gives 

that for all t > 0, 

Pr(Gt > x) < e~ r ' (p)x , X e <P (11) 
where r*(P) as defined in Assumption 1 is the positive root of the cumulant generating function of Z t . 

Lemma 1. Given infinite battery capacity, the battery level at the typical harvester satisfies 

S t <G t + G' t . (12) 

The proof of Lemma 1 is provided in Appendix A. Using (11) and Lemma 1, the threshold-crossing 
probability for the battery level can be shown to be bounded as follows. 

Lemma 2. Given infinite battery capacity and A e < P, the distribution of the battery level St satisfies 

Pv(S t >x)< 2e- 1 * r *( p K x - 2P \ (13) 

The proof of Lemma 2 is given in Appendix B. Define the energy-overshoot function D t : K, + — > M> + 
as the expected amount of energy stored in the typical harvester in excess of a threshold x > in the t-th 
slot: 

POD 

D t (x)= / (y-x)f s (y,t)dy (14) 

J X 

where f s (y,t) represents the probability density function of St. The function D t (x) can be bounded as 
shown in Lemma 3, which is proved in Appendix C using Lemma 2. 



> x 
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Lemma 3. Given infinite battery capacity and a threshold x > 2P, the energy-overshoot function satisfies 

A(x) < ^e-^ (P)( *- 2P) , Vt>0. 
Using Lemma 3, the main result of this section is readily obtained as shown below. 
Theorem 1. Given infinite battery capacity, the transmission probability is 



P= \ 



1, A e > P 

% A e <P. 



(15) 



Proof: First, consider the case of A e > P. Replacing the indicator function in (8) with one yields a 
lower bound on St, namely that St > Y2m=i ^m- As a result, p given in (4) can be lower bounded as 



1 n 

p > lim y E 

n^co n — ' 
t=l 



\m=l 



1 n (\ 1 P\ 

= lim -VPr - V(Z m -A e )>P-A e + - . (16) 

Using E[Z n ] = A e and applying the weak law of large numbers [25], for given r > and 5 > 0, there 
exists k > such that for all t > k, 

Pr^5^(Z m -Ae)>-rj >1-S. (17) 

Since A e > P , it follows that given 5 > 0, there exist r > and k! > such that r < (A e — P — j) for 
all t > k! . Using this fact and substituting (17) into (16), 



p > lim 1 V Pr f I V (Z m - A e ) > t J 

t=max(fe,fc') \ m=l / 

1 ™ 
> lim - V (1 - 5) 
n— >oo n z — ' 

t=max(fe,fc') 

[n-max(A;,fc / )](l-5) 



= lim 

= 1-5. (18) 

As S is arbitrary and p < 1, the result in the theorem statement for the case of A e > P follows from (18) 
by letting 5 — > 0. 

Next, consider the case of A e < P. The expected total amounts of harvested and transmitted energy by 

the i-th slot differ by the battery level in the t-th slot, namely 

t t 

E[Z m ] =PJ2 ^(Sm-i > P)] + E[5 t ]. (19) 

m=l m=l 
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Since E[S t ] > 0, 

1 * 

A e >Plim - VE[/(5 m _i >P)] 

t^>oo t L — ' 

m=l 

= Pp (20) 
where (20) substitutes (4). It follows that 

P<^, \e<P (21) 

Note that E[St] < x + D t (x) with x > 0. Using Lemma 3, for given 5 > 0, there exists x > such that 
E[Sy < x + 5. Combining this fact and (19) yields 

A e = P lim J V E[J(5 t _i > P)} + lim 

m=l 

< P lim I y E[I(S t -i > P)\ + lim ^ 

m=l 

= P P . 

As a result, 

p>^, X e <P. (22) 

The desired result for the case of A e < P is proved by combining (21) and (22), completing the proof. 

Last, the desired result for the boundary case of A e = P is proved by using the results proved above 
for A e 7^ P and letting P — > A e from either the right or the left, completing the proof. ■ 

A remark is offered as follows. According to Theorem 1, if P < A e , each transmitter transmits 
continuously in the steady state and is free of interruption caused by energy shortage. However, continuous 
transmissions are at the cost that the fraction of harvested energy at the rate of (A e — P) is never used 
for transmission and hence wasted. Next, if P > A e , there exists nonzero probability that the battery level 
of a transmitter is below P and hence transmission can be interrupted. Nevertheless, all harvested energy 
will be eventually used for transmission. 

B. Finite Battery Capacity 

The dynamics of the battery level at the typical harvester are characterized in (3). Let Dt(x) denote the 
energy-overshoot function for the case of finite-battery capacity, which is defined similarly as D t (x) in the 
preceding section. Given finite battery capacity, battery overflow can occur such that the battery saturates 
and arriving energy has to be discarded. The expected amount of discarded energy is measured by Dt{B) 
and can be bounded by the same upper bound on D t (x) (see Lemma 3) as shown below. 
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Lemma 4. Given finite battery capacity B > 2P and A e < P, the energy-overshoot function satisfies 

4 



D t (B) < 



r*(P) 



e -ir'(P)(B-2P)_ 



(23) 



Lemma 4 is proved in Appendix D. Using Lemma 4, bounds on the transmission probability are obtained 
as follows. 



Proposition 1. Given finite battery capacity B > IP, the transmission probability satisfies 



A e r*(P) 



e 2 



\r*{P){B-2P) 



<P< 



if \ e < P and 



if A e > P. 



max 

b>0 { A e + x 



A e r*(A e + x) 



e 2 



±r*{\ e +x)(B-2\ e -2x) 



< p < 1 



(24) 



(25) 



Proof: Consider the case of A e < P. By accounting for the discarded energy due to battery overflow, 
the expected total amounts of transmitted and harvested energy by the £-th slot is related by modifying 
(19) as 

t t t 

E[Zm] =PJ2 ^(Sm-l > P)} + E[S t ] + £ D t (B) 
m=l m=l m=l 

where the last term gives the expected amount of discarded energy. Applying Lemma 4 and St < B yields 

t 



A e < P lim - V E[J(S m _! > P)] + lim § + 



< Pp + 



m=l 

4 



t^oo t r*{P) 



e -\r*(P)(B-2P) 



r*{P) 



e -\r'{P){B-2P) 



and the first inequality in (24) follows. The second inequality is proved using Theorem 1 and the fact that 
limiting the battery capacity reduces p. 

Next, consider the case of A e > P. Let p'(x) denote the transmission probability for the virtual scenario 
where all transmissions use the power (P + x) with x > A e — P. It can be proved similarly as the first 
inequality in (24) that 



p'(x) > 



P + x 



1 



_ e -\r*(P+x)(B~2P~2x) 



X > A e - P. 



(26) 



A e r*(P + x) 

Note that p'{x) also gives the transmission probability for a virtual transmission strategy that removes x 
unit of energy from the battery of a harvester following every instance of transmission. Since removing 
energy from the battery reduces the transmission probability, p > p'{x) holds. Combining this inequality 
and (26) yields 



P> 



K 



P + x 



1 - 



\ e r*{P + x) 



e -\r*(P+x)(B-2P-2x) 



V X > A e - P. 
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It follows that 

A 

p > max 



I 4 c -h_r'(P+x)(B-2P-2x) 

X e r*(P + x) 



(27) 



x>\ e ~P P + X 

Manipulating the maximization variable x in (27) gives the first inequality in (25). The second inequality 
is obtained based on the same argument as the second inequality in (25), completing the proof. ■ 
Two remarks are offered as follows. 

1) For a sanity check, it can be observed from Proposition 1 that as B — > oo, p converges to its 
counterpart for the case of infinite battery capacity as stated in Theorem 1. 

2) It can be obtained from (24) that 

As observed from (28), the degradation of p due to finite battery capacity decreases exponentially 
with B. Hence the effect of finite battery capacity on p is expected to diminish rapidly as B increases 
if P is larger than A e , which agrees with simulation results (see Fig. 7). 
In general, exact analysis of p for the case of finite battery capacity is challenging except for some 
special cases, two of which are discussed as follows. 

1) Special case: bounded energy arrivals: Consider the case that Zt has bounded support and Zt € 
[0, 2 max ] . 

Proposition 2. Consider bounded energy arrivals. If z max < P and the battery capacity B > 2P, the 
probability for battery-overflow is zero and the transmission probability is 

P = y (29) 

where A e < P. 

Proof: By expanding (3), 

S t = 



mm(St-i + Z t -P,B), S t -i > P 

(30) 

mm(S t -i + Z t ,B), S t -i < P. 



Given B > 2P and Z t < P, it follows from (30) that S t -i < 2P ensures S t < 2P. Since S = 0, this 
result implies zero battery-overflow probability. Consequently, the transmission probability is identical to 
that for the case of infinite battery capacity in Theorem 1, proving the desired result in (29). ■ 
2) Special case: discrete energy arrivals: Assume that Z t is a discrete random variable and takes 
on values from {0, 1,2, • • • }, and that P is a positive integer. Under these assumptions, the distribution 
of battery levels can be analyzed using Markov-chain theory. Since St is independent with {S n } t r ^ 1 
given St-i, {St} satisfies the Markov property and is hence a Markov chain. Given battery capacity 
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Fig. 2. A Markov chain modeling the battery level at the typical transmitter for the special case of binary energy arrivals and 
P = 3. 



B, the Markov chain {St} has the state space {0, l,--- ,B}. Define the transition probability p rnn as 
Pmn = Pr (St = n I St-i = m). It can be obtained that 

Pr (Z t = n — m) , 
Pr (Z t = n - m + P) , 



Pr, 



n — m) 



m < P,n < B 
m > P,n < B 

m < P,n = B 



k=B—m 

oo 

Pr (Z t = n — m + P) , m > P,n = B. 



(31) 



k=B— m 

Let 7r m denote the steady-state probability of state m of the Markov chain {St}. Moreover, let P represent 
the transition-probability matrix with the (m, n)-th element given by p mn and 7r the steady-state-probability 
row vector with the m-th element given by ir m . Applying Perron-Frobenius theorem [25], 



7rP = 7T. 



(32) 



The stationary probabilities {7r m } can be computed by solving (32) under the constraint J2 m=0 vr m = 1. 
Given 7r, the transmission probability is obtained as p = Y^m=p n rn- 

The stationary probabilities {-ir m } can be derived in closed form for the simple case of binary energy 
arrivals, namely that Z t G {0, 1}. Then energy- arrival rate is A e = E[Z t ] = Pr(Z t = 1). The corresponding 
transition probabilities are modified from (31) as 

1 — A e , (m = P, n = 0) or (m < P, n = m) 



Pr, 



(33) 



(m = P, n = 1) or (m < P, n = m + 1) 
0, otherwise. 

Based on above transition probabilities, the Markov chain {St} is illustrated in Fig. 2. The stationary 
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distribution {ir m } satisfies the following equations obtained from (32) and (33): 

7T = 7Tp(l - A e ) + 7T (1 - A e ) 
7Tl = 1TpX e + 7Tl(l — A e ) + 7ToA e 

(34) 

7T m = VT m _iA e + 7T m (l - A e ), 771 = 2, ■ ■ ■ , P - 1 

TTp = 7Tp_iA e . 

Solving the equations in (34) and J2m=o n rn = I gives the following proposition. 

Proposition 3. Consider the case that energy arrivals are binary (Z t G {0, 1}) and the transmission power 

P is a positive integer no larger than the battery capacity B. The distribution of the battery level St is 

1 - A e A e 
vr = — p— , vtp = — 

7r m = p m = l,--- ,P- 1. 

The probability for battery overflow is zero and the transmission probability is p = irp = \ e /P. 

IV. Network Throughput 

In this section, using results on transmission probability as derived in the preceding section, the network 
throughput is maximized by optimizing transmission power under the outage -probability constraint and 
assuming infinite battery capacity. It is straightforward though tedious to extend the network-throughput 
analysis to the case of finite battery capacity using related results from the last section. 3 Such an extension 
provides few new insights and thus is omitted. 

A. Maximum Network Throughput 

To characterize the network throughput, transmission power P is related to the active transmitter density 
Xt under the outage-probability constraint. To this end, we define a parameter p e , called the nominal node 
density, as the density of a homogeneous PPP A(p £ ) such that 

Prl £ |T|- a >l) = e (35) 
VreA(M e ) / 

where the summation can be interpreted as the interference power measured at a receiver located at the 
origin from unit-power interferers distributed as A(p e ). Note that p e depends only on e and the distribution 



'Essentially, the network-throughput analysis for the case of finite battery capacity can be modified from that for the case of 
infinite battery capacity by replacing the transmission probability p with cp, where the factor c G (0, 1] represents the degradation 
due to limiting the battery capacity and can be quantified using results from Section III-B. 
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Fig. 3. The relation between the nominal node density [i a and the maximum outage probability e 

of a PPP and is independent with other network parameters. Moreover, fj, e is a strictly-monotone-increasing 
function of e due to the fact that denser interferers result in larger outage probability for a link. The 
expression of fi e has no closed form and its value can be computed by simulation (see e.g., [28]). The 
relation between fx e and e is shown in Fig. 3. The following lemma results from Mapping Theorem [26, 
pl8]. 

Lemma 5. The process aA(fj, t ) with a > is a homogeneous PPP with density /i e /a 2 . 

Define the admissible set F as all combinations of (At, P) that satisfy the outage-probability constraint: 

F = {(A*, P) g E+ x R+ | P out (A t , P) < e}. (36) 

A combination (At, P) is admissible if it belongs to F. To derive F, applying Lemma 5 allows P ut in 
(5) to be rewritten as 

Combining (35), (36) and (37) leads to the following lemma. 

Lemma 6. The admissible set F is given as 

-F^{(A t ,P)GlR + x]R+|A t <C(P)} (38) 
where C(P) represents the network interference temperature and is given as 

C(P) = K [I - p) ■ (39) 
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Two remarks are in order. 

1) The network interference temperature ((P) specifies the maximum density of interferers a link 
can tolerate without violating the outage-probability constraint. This quantity is analogous to the 
interference temperature in cognitive-radio systems that measures the maximum amount of addi- 
tional interference for a certain frequency band without significantly degrading the reliability of 
communications therein [29]. It can be observed from (39) that ((P) increases with growing P and 
reaches its maximum as P — > oo: 

lim C(P) = (40) 

2) The admissible set F is illustrated as the shaded region in Fig. 4. Increasing fi t (relaxing the outage- 
probability constraint) enlarges T and vice versa. The boundary of T corresponds to \ t = ((P). 

3) The condition X t < ((P) for the admissible set in (38) guarantees that P (equal to the received 
SNR) is no smaller than that is the minimum received SINR required for correct decoding. 

It is ready to derive the maximum network throughput R* and the optimal transmission power P*. 
The derivation can be intuitively explained using Fig. 4. It can be observed form the curve depicting the 
function Xt{P) that Xt(P) is fixed at Ao for all P < X e and a strictly-monotone-decreasing function for 
P > A e . Let P* correspond to the intersection between the curve (P, Xt(P)) and the admissible set T. 
Then under the outage-probability constraint, Xt(P) is maximized at P = P* , which corresponds to R* 
since it is proportionally to A^. Depending on if P* is larger or smaller than A e , P* and R* have different 
expressions as shown in Theorem 2 that states the main result of this section. 

Theorem 2. Given infinite battery capacity, the maximum network throughput R* and the optimal trans- 
mission power P* are given as follows. 

1) If Ao < C(A e ), 

R* = A log 2 (l + c?) (41) 



and P* is an arbitrary value in the range 
2) If Ao > C(A e ), 



1-6»(A /Me)2 



,A 



iT = ^log 2 (l + 0) (42) 



where P* > A e and VP* solves the following polynomial equation: 

x « _ 9 x a - 2 - 9 ' = 0. (43) 

Proof: First, consider the case of Ao < C(A e )- The throughput expression in (7) implies that 

R* < A o log(l + 0). (44) 
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Fig. 4. The admissible set T is plotted as the shaded region that contains all combinations of active transmitter density At and 
transmission power P that satisfy the outage-probability constraint. Given infinite battery capacity and fixed transmitter density 
Aq, At(-P) is a monotone-decreasing function of P as plotted with the thick line based on Theorem 1. Moreover, P* denotes the 
optimal transmission power that maximizes the network throughput. 



Define Pq such that Ao = ((Pq)- Using the definition of C in (39), 

i - e(A //. e )> 

It can be observed from (39) that ( is a strictly-monotone-increasing function. As a result, since Ao < C(^e) 
and Ao = C(^b)> Pq < and thus the set [Po, A e ] (the range of P* in the theorem statement) is nonempty. 
Consider an arbitrary value p € [i-foi^e]- Since Ao < C{p) from the monotonicity of 5, (Xt,P) = (Ao,p) 
is admissible according to Lemma 6. Furthermore, (X t ,P) = (\o,p) is feasible as p(p) = 1 based on 
Theorem 1 and p < A e . It follows that P = p maximizes R by achieving the equality in (44). This proves 
the desired result for the case of Ao < C(-^e)- 

Next, consider the other case of Ao > C(^e)- Using the throughput expression in (7) and Lemma 6, the 
problem of maximizing the network throughput is equivalent to 

maximize p{P) 

(45) 

subject to A p(P) < C(P). 
The inequality P* > A e in the theorem statement can be proved by contradiction as follows. Assume 
that P* < A e . This assumption results in p(P*) = 1 by applying Theorem 1. Moreover, < C(^e) 

is obtained using the aforementioned monotonicity of C> and hence Ao > C(P*) given that Ao > C(^e)- 
Combining Ao > C(P*) an(1 p(P*) = 1 shows that the assumption of P* < A e violates the constraint 
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in (45), proving that P* > A e . Then applying Theorem 1 yields the desired result in (42). Last, since 
p(P) and C(P) are strictly-monotone-decreasing and strictly-monotone-increasing functions, respectively, 
the solution P* for the problem in (45) must satisfy Xqp(P*) = ((P*) or equivalently \[P* solves the 
polynomial equation in (43), completing the proof. ■ 
Two remarks are in order. 

1) For the case of Ao < C(^e)> tne network is relatively sparse and A e is sufficiently large such that it 
is optimal as well as feasible for all transmitters to transmit with probability one, resulting in the 
network throughput in (41). For the case of Ao > C(^e)> the network is relatively dense and not all 
transmitters can transmit simultaneously under the outage -probability constraint even if they have 
sufficient energy for transmission. This corresponds to the optimal transmission probability smaller 
than one and the network throughput in (42). 

2) With a > 2 and the last coefficient at the right-hand side being negative, the polynomial equation 
in (43) has at least one strictly positive solution that gives P* for the case of Ao > ({X e ). For the 
special case of a = 4, the polynomial equation in (43) is quadratic and solving it gives P* in closed 
form as shown below. 

Corollary 1. Given infinite battery capacity, Aq > C(^ e ) and a = 4, the optimal transmission power P* is 



P* = — * 2 ^ J ■ (46) 

B. Maximum Network Throughput: Extreme Cases 

Consider a network with a high energy-arrival rate (A e — > oo). The maximum network throughput R* 
can be upper bounded as 

R* < ^log 2 (l + 0) (47) 

9 a 

since the outage-probability constraint requires that (see Lemma 6) 

A* < C(P) 

< h (48) 

9 a 

where (48) follows from (40) and that ( is a monotone-increasing function. Combining the two upper 
bounds on R* in (44) and (47) gives 

R* < min ( Ao, ^W 2 (l + (49) 



For a high energy-arrival rate, the equality in (49) is achieved as shown below. 
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Proposition 4. Given infinite battery capacity, as the energy-arrival rate A e — > oo, the maximum network 
throughput converges as 

lim R*(X e ) = min( A , ^ J log 2 (l + 0). (50) 



Proof: First, consider the case of Ao < Set P = A e — 5 with 5 > 0. This results in p = 1 

according to Theorem 1. Consequently, Ao = \t and hence Ao < // e a using (48). Combining this 
inequality and (40) yields that At < ((P) as P — > oo along with A e — > oo. It follows that as A e — > oo, the 
combination (X t , P) = (Ao, A e — 5) is admissible according to Lemma 6. This proves the equality in (49) 
for the current case. 

2 

Next, consider the case of Ao > n e Given this strict inequality, there exists 5 > such that 

Ao > ^ - S. (51) 



6-. 



Set P as 



(52) 



p _ A e Ap 

Combing (51) and (52) gives A e < P. Consequently, applying Theorem 1 gives 

fi e 6~l -5 

and hence 

Xt = 4 " S - (53) 

9 a 

As P — > oo along with A e — > oo, it follows from (40) and (53) that there exists 5 > such that At < C(-P)- 
As a result, by applying Lemma 6, the combination of (Xt,P) as specified in (52) and(53) is admissible. 
This leads to 

lim R*(X e ) = (^-5) log 2 (l + 0). (54) 

Letting 5—^0 proves the equality in (50) for the current case, completing the proof. ■ 
Two remarks are in order. 

1) Given a high energy arrival rate and infinite battery capacity, in the steady state, transmitters always 
have sufficient energy for transmission. Therefore, the expression in (50) also specifies the maximum 
network throughput of a MANET with reliable power supplies instead of energy harvesting. 

2 

2) For the case of Ao < °, the active transmitter density is below the network-interference 
temperature even though all transmitters transmit with probability one. Therefore, there is margin 
for further increasing active transmitter density without violating the outage-probability constraint. 
For this reason, the network-throughput limit in (50) for the current case is proportional to the 
transmitter density. However, for the case of Aq > Me/0"> active transmitter density reaches the 
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network-interference temperature and cannot be further increased. Consequently, the corresponding 
network-throughput limit in (50) is independent with the transmitter density. 
Consider a sparse network (Ao — > 0). It is optimal for each transmitter to transmit with probability 
one or equivalently set the transmission power as P* < A e (see Theorem 1). The corresponding network 
throughput is R* = Aolog(l + 6) that diminishes as Ao — > 0. Next, for a dense network (Ao — > oo), the 
maximum network throughput is specified as follow. 

Proposition 5. Given infinite battery capacity, as the transmitter density Ao — > oo, the maximum network 
throughput converges as 

lim iT(Ao) = ^log 2 (l + 9). (55) 
Proof: Set P as the following function of Ao: 

P^) = ^(\-^A~ l (56) 



H e \Q logAo 

that is shown shortly to achieve the limit of R* in (55). Given (56), there exists t\ > such that P(\o) > A e 
for all Ao > r. Therefore, it follows from (7) and Theorem 1 that 

lim R(X ,P(X )) = lim -^Ll O g 2 (l + 0) 

lim ^(I-^—V l og2 (l + 0) (57) 



A ->-oo \9 logAo 

= g-log 2 (l + 0). (58) 

Combining (47) and (58) shows that the maximum network throughput has the limit in (55) as Ao increases. 

The remaining proof verifies that P(\q) and the corresponding X t are admissible as Ao — > oo. It follows 
from (57) that for all Ao > t\, Xt is a function of Ao and given as 

a ' <a °»="-0-i^)"- <59) 

Substituting (56) into (39) yields 

C(P(A )) = ^(\-^(\- AA ■ (60) 



6 A A e \9 logA 0/ 

By comparing (59) and (60), there exists T2 > such that At(Ao) < ((P(Xq)) for all Ao > t 2 . This proves 
the admissibility of P(Xq) in (56) and X t (Xo) in (59) as Ao — > oo, completing the proof. ■ 
A remark is offered as follows. The rate of total energy harvested per unit area is A^P* = AoA e as 
Ao — > oo. The linear growth of the rate with increasing Ao is due to that the harvester density is equal to 
Aq. However, more aggressive energy harvesting in a dense network does not continuously increase the 
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network throughput that saturates at high transmitter power as the network becomes interference limited 
(see Proposition 5). This issue may be resolved by using an alternative multiple-access protocol such as 
frequency-hopping multiple access that reduces the density of simultaneous transmitters. 

V. Numerical Results 

The nominal node density fi e is fixed as 0.05 for all numerical results, corresponding to the maximum 
outage probability e « 0.015. The relation between fj, e and e is obtained by simulation based on the 
following procedure (see e.g., [28]). The summation over the PPP A(/i e ) in (35) is approximated by the 
signal power measured at the origin due to unit-power transmissions by transmitters uniformly distributed 
in a disk. The number of transmitters follows the Poisson distribution with mean 200 and the disk radius 
is adjusted such that the expected transmitter density is equal to /v Based on this setup, the values of 
(e, |U e ) are computed using the Monte Carlo method that yields the plot in Fig. 3. In addition, all numerical 
results are based on the SINR threshold 8 = 3 and the path-loss exponent a = 3. 

Infinite battery capacity is assumed for the numerical results presented in Fig. 5 and Fig. 6. In Fig. 5, 
the maximum network throughput R* computed using Theorem 2 is plotted against the increasing energy- 
arrival rate A e for the transmitter density Ao = {0.02,0.05,0.5}. It can be observed from Fig. 5 that R* 
grows as A e increases and saturates for large A e . The limits agree with those computed using Proposition 4, 
namely 0.04 bit/s/Hz/unit-area for A = 0.02 and 0.048 bit/s/Hz/unit-area for A = {0.05, 0.5}. In addition, 
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Fig. 6. Maximum network throughput versus transmitter density for optimal transmission power, infinite battery capacity, and 
the energy-arrival rate A e = {0.5, 1,5}. 

Fig. 5 shows that in a denser network (i.e., A e = 0.5), R* reaches its limit more rapidly as A e increases. 

Fig. 6 shows the curves of R* versus Ao for A e = {0.5, 1, 5}, which are obtained using Theorem 2. As Ao 
increases and regardless of the value of A e , R* is observed to converge to the limit 0.048 bit/s/Hz/unit-area 
predicted by Proposition 5. Moreover, it is observed from Fig. 6 that larger A e results in faster convergence 
of R* to its limit as Ao increases. 

Next, consider the case of finite battery capacity. To investigate the corresponding transmission probabil- 
ity and maximum network throughput, the distribution of the energy arrival process {Z t } is particularized 
as follows. Let {Vt} denote an i.i.d. sequence of random variables following the chi-squared distribution 
with d G {1,2, •••} degrees of freedom (DoF) and mean equal to d\ e . Let {Zt} = {^Vt} and hence 
Zt has mean A e and variance 2X^/d. The chosen distribution of Z t allows its variance (randomness) to 
be controlled by varying d while the mean of Z t is fixed. Note that Z t converges to a constant A e in 
probability as d — > oo by the law of large numbers. 

In Fig. 7, transmission probability p obtained by simulation is plotted against increasing transmission 
power P for finite battery capacity B = {1.5P, 2P, 4P, 10P}, d = 4, A e = 2, and A = 0.02. It is found 
that B = 10P is sufficiently large such that the values of p closely match those for the case of infinite 
battery capacity as computed using Theorem 1. As observed from Fig. 7, finite battery capacity degrades 
p significantly only when P and hence B are relatively small; as P and B increase, p rapidly approaches 
the counterpart for the case of infinite battery capacity (or that for B = 10P). 
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Fig. 7. Transmission probability versus transmission power for finite battery capacity B — {1.5P, 2P, 4P, 10P}, the DoF of the 
energy-arrival process d — 4, the energy-arrival rate A e = 2, and the transmitter density Ao = 0.02. 



Fig. 8 displays the curves of p versus P obtained by simulation for B = 1.5P, d = {2,4,8,16}, 
A e = 2 and Ao = 0.02. For comparison, the curve for the case of infinite battery capacity is also plotted. 
As observed from Fig. 8, reducing the randomness of the energy arrival process by increasing d leads to 
smaller battery-overflow probability and hence higher p. The effect of d on p diminishes as P (and hence 
battery capacity) increases and p converges to its counterpart for the case of infinite battery capacity. 

Last, we investigate the effect of the DoF of the energy arrival process on the network throughput. To 
this end, Fig. 9 shows the curves of R* versus A e obtained by simulation for B = 1.5P, d = {2, 4, 8, 16}, 
and Ao = 0.02. The curve for the case of infinite battery capacity is also plotted for comparison. It can 
be observed from Fig. 9 that finite battery capacity causes significant throughput loss especially for large 
A e . Such loss is smaller for larger d because of less randomness in harvested energy and hence smaller 
battery-overflow probability. 

VI. Conclusion 

The energy dynamics in a mobile ad hoc network have been characterized in terms of transmission 
probability. Assuming infinite battery capacity, it has been found that the transmission probability is equal to 
one when the energy-arrival rate exceeds transmission power or otherwise is equal to their ratio. Moreover, 
for the case of finite battery capacity, bounds on the transmission probability have been obtained and exact 
expressions have been derived for the special cases of bounded or discrete energy arrivals. The results on 
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Fig. 8. Transmission probability versus transmission power for both the cases of finite (B — 1.5P) and infinite battery capacity. 
The DoF of the energy-arrival process is d = {2, 4, 8, 16}, the energy- arrival rate X e — 2, and the transmitter density Ao = 0.02. 



transmission probability have been applied to derive the maximum network spatial throughput for a given 
energy-arrival rate and optimized transmission power. It has been shown that it is optimal and feasible for 
all transmitters to transmit with probability one if the transmitter density is below a threshold that depends 
on the energy- arrival rate; otherwise, each transmitter should transmit with probability smaller than one. 

There are several potential directions for extending this work. Coexisting wireless networks may harvest 
electromagnetic (EMR) energy from each others' transmissions. Modeling and designing coexisting net- 
works with EMR energy harvesting give rise to many new research issues ranging from algorithm design 
to throughput analysis. The current work focuses on ad hoc networks with random access and can be 
extended to other types of networks such as cellular networks or other medium-access-control protocols 
such as carrier- sensing multiple access. Last, it is interesting to investigate the effects of bursty data arrivals 
and more sophisticated power control on the throughput of wireless networks with energy harvesting. 
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Appendix 

A. Proof of Lemma 1 

If t 6 T, it follows from (10) that St = G' t and hence the inequality in the lemma statement holds since 
Gt > 0. Next, consider the case of t £ T and let to 6 T denote the time instant closest to but smaller 
than t. Then St can be upper bounded as 

S t <m^(0,S to+ ^l to+i Z t ) 

<^ + max(0,^ o+i Z,) 

<Sfc+max(0,Zt,Z t + Z t _i,--- ,V* Z t ) . (61) 

It can be obtained from (9) that 

G t = max (o, Z t , Z t + Z^, ■ ■ ■ , V* , ^ Z u G to +V* Z t ) . (62) 

Since G to > and G' ta = St from (10), combining (61) and (62) proves the inequality in the lemma 
statement for the case of t £ T, completing the proof. ■ 
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B. Proof of Lemma 2 

Using Lemma 1 and for < a < x, 

Pr(Si > x) < Px{G t + G' t >x) 

= Pi{G t > x - G' t | G' t > a) Pr (G' t > a) + Pr(G t > x - G\ \ G' t < a) Pr (G' t < a) 

< Pr (G' t >a)+ Pr(G t > x - a). (63) 

Let to £ T specify the slot such that G' t = St . Since <St _i < P based on the definition of T, St = 
St -i + Zt„ using (8). It follows from this equality and (63) that 

Pr{S t >x)<Pt (S to -i + Z to >a) + Pr(G t > x - a) 
< Pr (Z to >a-P) + Pv(G t > x - a) 

= Pr(Z to >a-2P)+Pr(G t >x-a). (64) 
By bounding the first term in (64) using Chernoff bound [25] and the second using (11), 



Pr(St > x) < min E 



e r Zl 



e -r(a-2P) +e -r'(P)(x-a) 



r>0 

< e -r-(P)(a-2P) + r «(P)(,- a) (65) 

where (65) results from setting r = r*(P). By choosing a such that the exponents of the two terms in 
(65) are equal, the desired result follows. ■ 

C. Proof of Lemma 3 

From the definition in (14), 

POD 

D t (x) = yfs(y,t)dy - xPr(S t > x) 

J X 



Pr(S t > y)dy (66) 

X 

poo 

< / 2e-^' (p ^- 2P Uy (67) 

J X 

where (66) is obtained using integration by parts. The desired result follows from (67). ■ 
D. Proof of Lemma 4 

Given finite battery capacity B, the inequality St < G' t + Gt can be proved using induction as follows. 
This inequality holds for t = since S = G = and G' Q = P. Assume that S t < G' t + G t . If t + 1 e T, 
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't+i 



< < 



since G' t+1 = St+i from the definition of {G' t } in (10), 

St+i < St+i + min(max(G ( + Z t+1 , 0), B — P) 

= G' t+1 + G t+1 . (68) 
If t + 1 ^ T, based on the evolution of {St} in (3), 

mm{S t + Z t+1 ,B), S t < P 
mm{S t + Z t+1 ,B), S t > P 
P, S t <P 
S t + min(max(G t + Z t+1 , 0) , B) , S t > P 
< G' t+1 + G t +i (69) 

where (69) uses G' t+1 > P and the definition of {Gt} in (9). Combining (68) and (69) proves that 
S t +i < G' t+1 + G t +i if St < G[ + G t . It follows that S t < G' t + G t for all t > 0. Using this inequality 
and following the same procedures as for proving Lemma 2 and Lemma 3, it can be shown that (13) also 
holds for the case of finite battery capacity and 

A(x) < -±-e-\r*W*-*F>, Vt>l. (70) 

The desired result follows by setting x = B. ■ 



